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In this paper, for the first time a method is proposed to compute electromagnetic 
effects in hadronic processes using lattice simulations. The method can be applied, 
for example, to the leptonic and semileptonic decays of light or heavy pseudoscalar 
mesons. For these quantities the presence of infrared divergences in intermediate 
stages of the calculation makes the procedure much more complicated than is the 
case for the hadronic spectrum, for which calculations already exist. In order to 
compute the physical widths, diagrams with virtual photons must be combined with 
those corresponding to the emission of real photons. Only in this way do the infrared 
divergences cancel as first understood by Bloch and Nordsieck in 1937. We present a 
detailed analysis of the method for the leptonic decays of a pseudoscalar meson. The 
implementation of our method, although challenging, is within reach of the present 
lattice technology. 
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I. INTRODUCTION 


Precision flavour physics is a particularly powerful tool for exploring the limits of the 
Standard Model (SM) of particle physics and in searching for inconsistencies which would 
signal the existence of new physics. An important component of this endeavour is the over¬ 
determination of the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix from a 
wide range of weak processes. The precision in extracting CKM matrix elements is generally 
limited by our ability to quantify hadronic effects and the main goal of large-scale simula¬ 
tions using the lattice formulation of QCD is the ab-initio evaluation of the non-perturbative 
QCD effects in physical processes. The recent, very impressive, improvement in lattice com¬ 
putations has led to a precision approaching 0(1%) for a number of quantities (see e.g. 
Ref. [lj and references therein) and therefore in order to make further progress electromag¬ 
netic effects (and other isospin-breaking contributions) have to be considered. The question 
of how to include electromagnetic effects in the hadron spectrum and in the determination 
of quark masses in ab-initio lattice calculations was addressed for the first time in 0. Much 
theoretical and algorithmic progress has been made following this pioneering work, particu¬ 
larly in recent years, leading to remarkably accurate determinations of the charged-neutral 
mass splittings of light pseudoscalar mesons and light baryons (see Refs. ( 3 - 8 ] for recent 


papers on the subject and Refs. 


fl 


IP] for reviews of these results and a discussion of the 


different approaches used to perform QED+QCD lattice calculations of the spectrum). 

In the computation of the hadron spectrum there is a very significant simplification in 
that there are no infrared divergences. In this paper we propose a strategy to include elec¬ 
tromagnetic effects in processes for which infrared divergences are present but which cancel 
in the standard way between diagrams containing different numbers of real and virtual 
photons ll|. The presence of infrared divergences in intermediate steps of the calculation 
requires the development of new methods. Indeed, in order to cancel the infrared divergences 
and obtain results for physical quantities, radiative corrections from virtual and real pho¬ 
tons must be combined. We stress that it is not sufficient simply to add the electromagnetic 
interaction to the quark action because amplitudes with different numbers of real photons 
must be evaluated separately, before being combined in the inclusive rate for a given process. 
In this paper for the first time we introduce and discuss a strategy to compute electromag¬ 
netic radiative corrections to leptonic decays of pseudoscalar mesons which can then be used 








3 


to determine the corresponding CKM matrix elements. Although we present the explicit 
discussion for this specific set of processes, the method is more general and can readily be 
extended to generic processes including, for example, to semileptonic decays. 

We now focus on the leptonic decay of the charged pseudoscalar meson P + . Let r 0 be the 
partial width for the decay P + —> £ + vi where the charged lepton l is an electron or a muon 
(or possibly a r) and U£ is the corresponding neutrino. The subscript 0 indicates that there 
are no photons in the final state. In the absence of electromagnetism, the non-perturbative 
QCD effects are contained in a single number, the decay constant fp, defined by 

(0 | ?i7 m 7 5 <Z2 | P + (p)) — ip^fp , (1) 


where P + is composed of the valence quarks q\ and g 2 , and the axial current in (JT]) is 
composed of the corresponding quark fields. There have been very many lattice calculations 


of the decay constants /x,/ d ( s) and f B , 




EL 


some of which are approaching 0(1%) 


precision. As noted above, in order to determine the corresponding CKM matrix elements 
at this level of precision isospin breaking effects, including electromagnetic corrections, must 
be considered. It will become clear in the following, and has been stressed in jl2, lfj, that 
it is not possible to give a physical definition of the decay constant fp in the presence 
of electromagnetism, because of the contributions from diagrams in which the photon is 
emitted by the hadron and absorbed by the charged lepton. Thus the physical width is not 
just given in terms of the matrix element of the axial current and can only be obtained by 
a full calculation of the electromagnetic corrections at a given order. 

The calculation of electromagnetic effects leads to an immediate difficulty: T 0 contains 
infrared divergences and by itself is therefore unphysical. The well-known solution to this 
problem is to include the contributions from real photons. We therefore define Ti(AE) to 
be the partial width for the decay P + —y £ + v , £ r y where the energy of the photon in the 
rest frame of P + is integrated from 0 to A E. The sum T 0 + Ti(AE) is free from infrared 
divergences (although, of course, it does depend on the energy cut-off A E). We restrict the 
discussion to 0(a) corrections, where a is the electromagnetic fine-structure constant, and 
hence only consider a single photon. 

The previous paragraph reminds us that the determination of the CKM matrix elements 
V qiq2 at 0(a) (i.e. at 0(1%) or better) from leptonic decays requires the evaluation of 
amplitudes with a real photon. The main goal of this paper is to suggest how such a 
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calculation might be performed with non-perturbative accuracy. There are a number of 
technicalities which will be explained in the following sections, but here we present a general 
outline of the proposed method. We start with the experimental observable T(A E), the 
partial width for P + —> £ + ui( 7 ). The final state consists either of £ + U£ or of £ + v^ where 
the energy of the photon in the centre-of-mass frame is smaller than A E\ 


T(AE) = r 0 + T 1 (AE). 


( 2 ) 


In principle at least, Ti(AE) can be evaluated in lattice simulations by computing the 
amplitudes for a range of photon momenta and using the results to perform the integral 
over phase space. Such calculations would be very challenging. Since the computations are 
necessarily performed in finite volumes the available momenta are discrete, so that it would 
be necessary to choose the volumes appropriately and compute several correlation functions. 
We choose instead to make use of the fact that a very soft photon couples to a charged hadron 
as if to an elementary particle; it does not resolve the structure of the hadron. We therefore 
propose to choose A E to be sufficiently small that the pointlike approximation can be used 
to calculate T 1 (A E) in perturbation theory, treating P + as an elementary particle. On the 
other hand, A E must be sufficiently large that T(AE) can be measured experimental 
imagine setting A E = O(10-20MeV) which satisfies both requirements. From Refs, 
we learn that resolutions on the energy of the photon in the rest frame of the decaying 
particle of this order are experimentally accessible. In Appendix |B] we present a discussion, 
based on phenomenological analyses, of the uncertainties induced by treating the meson as 
elementary as a function of A E. 

It is necessary to ensure that the cancellation of infrared divergences occurs with good 
numerical precision leading to an accurate result for T(AE). Since To is to be calculated in 
a Monte-Carlo simulation and Ti(Aif) in perturbation theory using the pointlike approxi¬ 
mation, this requires an intermediate step. We propose to rewrite Eq. (J2]) in the form 



r(A-E) = lim (r„ - rg l ) + lim (rg* + r,(AB)). (3) 

V—too V —^oo 

where V is the volume of the lattice. r[f is an unphysical quantity; it is the perturbatively 
calculated amplitude at 0(a) for the decay P + —> £ + Uf with the P + treated as an elementary 
particle. In Tq’ the finite-volume sum over the momenta of the photon is performed over 
the full range. The contributions from small momenta to T 0 and Fq 1 are the same and thus 
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the infrared divergences cancel in the first term on the right-hand side of Eq. (j3]l . Moreover, 
the infrared divergences in r 0 and Tq* are both equal and opposite to that in Ti(AE). 
The infrared divergences therefore cancel separately in each of the two terms on the right- 
hand side of Eq. (J3J) and indeed we treat each of these terms separately. T pt + Ti(AE) is 
calculated in perturbation theory directly in infinite volume. The QCD effects in To are 
calculated stochastically in a lattice simulation and the virtual photon is included explicitly 
in the Feynman gauge. For each photon momentum this is combined with T pt and the 
difference is summed over the momenta and then the infinite-volume limit is taken. This 
completes the sketch of the proposed method, and in the remainder of this paper we explain 
the many technical issues which must be addressed. 

It will be helpful in the following to define AFo(T) in terms of the first term on the 
right-hand side of Eq. ([3]) : 

Ar„(L) = r„(L) - rs‘(L), ( 4 ) 

where we have made the dependence on the volume explicit, V = L 3 and L is the length of 
the lattice in any spacial direction (for simplicity we assume that this length is the same in 
all three directions). In analogy to Eq. ([2]) we also define the perturbative quantity 

r pt (A£) = r[? + ri(A£). (5) 

We note that, since the sum of all the terms in Eq. (j3J) is gauge invariant as is the 
perturbative rate r pt (AE), the combination Ar 0 (T) is also gauge invariant, although each 
of the two terms is not. 

The plan of this paper is as follows. In the next section we discuss the effective weak 
Hamiltonian and its renormalisation in the presence of electromagnetism. The structure 
of the calculation and the correlation functions which need to be calculated are presented 
in Sec. m The evaluation of the second term on the right-hand side of Eq. (l3lh r pt (AE), 
directly in infinite volume, is theoretically straightforward and we perform this calculation 
in Sec.lVl Sec. lVIl contains a detailed discussion of the regularisation and cancellation of 
infrared divergences in a finite volume. We put all the elements of the calculation together 
in Sec. lVIl! where we present a summary and the prospects for the implementation of the 
method in numerical simulations. There are two appendices. In AppendixlAl we discuss 
the matching of the bare lattice operators used in the calculation of correlation functions 
and those defined in the W -regularisation which is a natural scheme used in the definition 
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FIG. 1: Tree-level diagram for the process ud —>• l + vi (left-hand diagram). In the effective theory 
the interaction is replaced by a local four-fermion operator (right-hand diagram). 

of the Fermi constant Gf in the presence of electromagnetism. Finally in Appendix iBl we 
present some phenomenological estimates of the uncertainties due to the use of the point-like 
approximation for P + in the decay P + —y £ + uj. 

In the remainder of the paper, to be specific we choose P + = n + but the discussion 
generalizes trivially to other pseudoscalar mesons with the obvious changes of flavour labels. 
The method does not require P + to be a light psuedo-Goldstone Boson nor on the use of 
chiral perturbation theory. 


II. MATCHING THE EFFECTIVE LOCAL FOUR-QUARK OPERATOR(S) 

ONTO THE STANDARD MODEL 

At lowest order in electromagnetic (and strong) perturbation theory the process ud —y 
£ + ui proceeds by an s-channel W exchange, see the left-hand diagram in Fig.[Q Since the 
energy-momentum exchanges in this process are much smaller than Mw, it is standard 
practice to rewrite the amplitude in terms of a four-fermion local interaction: 

AG 

£w = ——^=- Vud { d LlfiU L ) ('A lY'£l) , (6) 

where the subscript L represents left , Vt = 2 75 and Gp is the Fermi constant. In 

performing lattice computations this replacement is necessary, since the lattice spacing a is 
much greater than 1 /Mw, where Mw is the mass of the lU-Boson. When including the 0(a) 
corrections, the ultra-violet contributions to the matrix element of the local operator are 
different to those in the Standard Model and in this section we discuss the matching factors 
which must be computed to determine the 0(a) corrections to the n + —y £ + ui decay from 
lattice computations of correlation functions containing the local operator in (TB1) . Since 
the pion decay width is written in terms of Gf, it is necessary to start by revisiting the 
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FIG. 2: Diagrams contributing to the 0(a ) corrections to muon decay; see Eq. (J7|. The curly line 
represents the photon. 

determination of the Fermi constant at 0(a). 

A. Determination of the Fermi constant, Gf 


Gf is conventionally taken from the measured value of the muon lifetime using the ex- 

fllrz] 
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27T \ 4 J 


( 7 ) 


leading to the value Gp = 1.16634 x 10 5 GeV 2 . (For an extension of Eq. (J7J) to 0(a 2 ) 
and the inclusion of higher powers of p = (m e /m^) 2 see Sec. 10.2 of 18| . The Particle Data 


18] quote the corresponding value of the Fermi constant to be Gp = 1.1663787(6) x 


Group 
10" 5 GeV" 2 .) 

Eq. ([7j can be viewed as the definition of Gf. When calculating the Standard Model 
corrections to the muon lifetime many of the contributions are absorbed into Gf and the 
remaining terms on the right-hand side of (jTJ) come from the diagrams in Fig. El Specifically 
in these diagrams the factor 1/k 2 in the Feynman-gauge photon propagator is replaced by 
1/k 2 x M^ V /(M( V — k 2 ), where k is the momentum in the propagator; this is called the 
W -regularisation of ultra-violet divergences. These diagrams are evaluated in the effective 
theory with the local four-fermion operator (G/Y m (1 — 7 5 )/i) (e 7 M (l — 7 5 )zy=); the two currents 
are represented by the filled black circles in Fig. [2] 

An explanation of the reasoning behind the introduction of the W-regularisation is given 


m 


19J. The Feynman-gauge photon propagator is rewritten as two terms: 

1 M 2 1 


1 

k 2 


k 2 -M 2 w 


M 2 , - k 2 k 2 


( 8 ) 


and the ultra-violet divergent contributions come from the first term and are absorbed in 
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FIG. 3: Photon-IF box diagrams contributing to the O(a) corrections to muon decay in the 
Standard Model. The curly line represents the photon. 

the definition of Gp. In addition, the Standard-Model 7 -W box diagram in Fig.[3]is ultra¬ 
violet convergent and is equal to the corresponding diagram in the effective theory (i.e. the 
third diagram in Fig. [2]) with the W-regnlarisation, up to negligible corrections of 0(q 2 /M^ v ), 
where q is the four-momentum of the electron and its neutrino. Other electroweak corrections 
not explicitly mentioned above are all absorbed into Gp. 


B. VF-regularisation and Weak Decays of Hadrons 


It is a particularly helpful feature that most of the terms which are absorbed into the 
definition of Gp are common to other processes, including the leptonic decays of pseudoscalar 
mesons { 20 , 21]. There are however, some short-distance contributions which do depend on 
the electric charges of the individual fields in the four-fermion operators and these lead to a 
correction factor of (1 + 2s l 0 g Mz.) to r 0 20]. This is a tiny correction (~ 0.06%), but one 
which nevertheless can readily be included explicitly. 

The conclusion of the above discussion is that the evaluation of the amplitude for the 
process tt + —y £ + u up to 0(a) can be performed in the effective theory with the effective 
Hamiltonian 


H cS = v* d ^1 + ^ log (g? 7 m (1 - 7 5 )w) (z77m(1 ~ 7 5 )^) , (9) 

and with the Feynman-gauge photon propagator in the W-regularisation. The value of Gf 
is obtained from the muon lifetime as discussed around Eq. ( 171 ) . 

Of course we are not able to implement the W-regularisation directly in present day lattice 
simulations in which the inverse lattice spacing is much smaller than M\p- The relation 
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between the operator in eq. (J9]) in the lattice and W regularisations can be computed in 
perturbation theory. Thus for example, with the Wilson action for both the gluons and 
fermions: 

0^" reg = (l + |-(2 log a 2 M^~ 15.539)) + £- (0.536 O h 2 we 

+1.6070^ are - 3.214 0\ are - 0.804 0£ are ) , (10) 

where 

Oi = (d 7 M (l - 7 b )u) (Pr +(1 - 7 b )£) 0 2 = (d 7 M (l + 7 5 )u) (Pr^(l - 7 5 )f) 

0 3 = (J(l - 7 5 )m) (P £ (l + 7 5 )£) 0 4 = (d(l + 7 5 )w) (P*(l + 7 5 )£) (11) 

0 5 = (Ja^(l + 7 5 )m) (z^ct M j,(1 + 7 5 )£). 

The superscript “bare” indicates that these are bare operators in the lattice theory and the 
presence of 5 operators on the right-hand side of Eq. (fTO is a consequence of the breaking of 
chiral symmetry in the Wilson theory. Using lattice actions with good chiral symmetry, such 
as domain wall fermions with a sufficiently large fifth dimension, only 0\ &ve would appear on 
the right-hand side of Eq. dTOlh The coefficients multiplying the operators depend of course 
on the lattice action being used. More details of the derivation of Eq. (fTOl) are presented in 
AppendixlAl Eq. (TTOT) is valid up to corrections of 0(a s (a ) a). 

Having formulated the problem of calculating To in terms of the evaluation of correlation 
functions involving the effective Hamiltonian in Eq. (J9J) we are now in a position to discuss 
the calculation of ATo(T), the first term on the right-hand side of the master formula Eq. (|3l) . 


III. STRUCTURE OF THE CALCULATION 

In this section we begin our explanation of how the calculations of the amplitudes for 
the processes n + —> £ + u and n + —> are to be performed. Before entering into the 

details however, we discuss more extensively the structure of the different terms appearing 
in Eq. ([3]) . 

Since we add and subtract the same perturbative quantity To*, we fold it convenient 
to choose this to be the virtual decay rate for a point-like pion computed in the W- 
regularisation. In this way we obtain the important advantage that the difference of the 
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i + 




FIG. 4: Correlation function used to calculate the amplitude for the leptonic decay of the pion in 
pure QCD. The two black filled circles represent the local current-current operator (07/A)i 

the circles are displaced for convenience. 

first two terms (AT 0 (L)) and the sum of the last two terms (T pt (AE)) on the r.h.s. of 
Eq. (J2D are separately ultraviolet and infrared finite. 

Let \J~Zi be the contribution to the decay amplitude from the electromagnetic wave- 
function renormalisation of the final state lepton (see the diagram in Fig.[5](d)). An impor¬ 
tant simplifying feature of this calculation is that Zi cancels in the difference To — r pt . This 
is because in any scheme and using the same value of the decay constant f n , the contribu¬ 
tion from the diagram in Fig.[5](d) computed non-perturbatively or perturbatively with the 
point-like approximation are the same. Thus we only need to calculate Z(_ directly in infinite 
volume and include it in the second term on the right-hand side of Eq. (I3lh As a result of 
this cancellation it is convenient to rewrite T 0 and Tf 1 in the form: 


p _ -ptree , pc 

1 0 — 1 n +b 


yd) 


and r pt = rff 66 + c pt + r, 


(d),pt 


( 12 ) 


where the superscript tree indicates the width in the absence of electromagnetic effects, (d) 
denotes the contribution from the leptonic wave function renormalisation and the index a 
represents the remaining contributions of 0(a) other than those proportional to Ze. In this 
notation the above discussion can be summarised by saying that Tg d) = Tg ,pt and that the 
calculation of AT 0 (L) at 0(a) reduces to that of computing Tq — Tg’ pt . 

Having eliminated the need to include the effects of the lepton’s wave-function renormal¬ 
isation from the evaluation of ATo(L), we need to make the corresponding modification in 
the factor(s) relating the lattice and W regularisations. This simply amounts to subtracting 
the term corresponding to the matching between the lattice to W regularisations of the lep¬ 
ton wave function renormalisation diagram. With the Wilson action (for both gluons and 
fermions) for example, the 0(a) contribution to this matching factor is 

Az w- r =g = X- (A - log a 2 MS, - 11.852 


(13) 
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Thus, with the Wilson action, we can avoid calculating the effects of the lepton’s wave- 
function renormalisation in AT 0 (T) by neglecting the diagram in Fig.^d) and the corre¬ 
sponding diagram with the point-like pion, and simply replacing Q^-reg in Eq. dlUli by 


£)W-reg 


(i + ^ log a 2 M^ - 8.863^ 0\ we + (0.536 O dare 

+1.607 0^ are - 3.214 0^ re - 0.804 0]? are ) . 


(14) 


Such matching factors depend, of course, on the lattice discretisation of QCD and we simply 
present the results for the Wilson action for illustration. 

Of course rQ d ' ),pt needs to be computed for the second term on the right-hand side of 
Eq. (ED- This is a straightforward perturbative calculation in inhnite-volume and gives 


~p(d),pt _ -ptree OL 

iQ _i ° 4 vT 


log 


m £ 

Wv 


- 2 log 


( m. 


\ m 


(15) 


where we use the VF-regularisation for the ultra-violet divergences and have introduced a 
mass m 1 for the photon in order to regulate the infrared divergences. The explicit expression 
for Tjf 6 is given in Eq. (I20l) below. Using the IF-regularisation we naturally work in the 
Feynman gauge, but note that with m 1 as the infrared regulator the result for Z\ is gen¬ 
erally gauge-dependent. For example, using dimensional regularisation for the ultraviolet 
divergences and m 7 as the infrared regulator leads to a gauge dependent result for this single 
diagram (gauge invariance is restored of course for r pt (AE)). 

In summary therefore, we need to compute the two quantities 


Ar 0 (L) = f“ - r“’ pt and r pt (A£) = r|, ree + r“’ pt + r[, d) ’ pt + t^ae) , (16) 

where fg corresponds to Tq using o^-i-eg j ns ^ eac j Q f ^w-re g ^ 0 |- e that Ar 0 (L) and r pt (AE) 
are separately infrared finite and the result of the calculation of these two quantities does 
not depend on the infrared cutoff. In particular, this means that the infrared cutoff can be 
chosen in two different ways for the two quantities. We have decided to give a mass to the 
photon in the perturbative calculation of r pt (AE), whereas for AFo(T) a possible convenient 
choice is to use the finite volume as the infrared regulator. This will be explained in more 
detail in Sec. ED 

In the following two sections we discuss the calculation of Ar 0 (T) and r pt (Ai7) respec¬ 
tively. 
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IV. CALCULATION OF Ar 0 (L) 

In this section we describe the calculation of the first term on the right-hand side of 
Eq. (Ql), AT 0 (L), at 0(a). We start however, by briefly recalling the calculation of r 0 at 
0(a°), i.e. without electromagnetism. 


A. Calculation of Tq at O(a 0 ) 


Without electromagnetic corrections we need to compute the correlation function 
sketched in Fig.[4l which is a completely standard calculation. Since the leptonic terms 
are factorized from the hadronic ones, the amplitude is simply given by 


G 

Uvta(Pvt)(M 0 ) a pv e p(p £ ) = -^j=V* d (Q\dYl b u\'K + (p^)) [u Ue (p Ue )^(l ~7 5 )ve(Pi)] 


^Gpf-n 


KdPn K,fc,b,.(i --y 6 )»/(?/)] 


(17) 


Here u, d in the matrix element represent the quark holds with the corresponding flavour 
quantum numbers and u Vf and the spinors of the leptons defined by the subscript. The 
hadronic matrix element, and hence the decay constant f n , are obtained in the standard 
way by computing the correlation function 


Co(t) = 53 <0 | ( rf (0, °bV «(0,0)) 0 f (£, -t) |0) ~ m °* A 0 , (18) 

x m ' K 

where 0f is an interpolating operator which can create the pion out of the vacuum, = 
(7r + (0)|<b(0, 0) | 0) and Aq = ( 0 | d'-f A 'y 5 u |7r + (0 )) 0 . We have chosen to place the weak current 
at the origin and to create the pion at negative time — t, where t and T — t are sufficiently 
large to suppress the contributions from heavier states and from the backward propagating 
pions (this latter condition may be convenient but is not necessary). The subscript or 
superscript 0 here denotes the fact that the calculation is performed at O(a 0 ), i.e. in the 
absence of electromagnetism. Zq is obtained from the two-point correlation function of two 
0 operators: 

Cf(i) = X) (19) 

x n 

For convenience we take 0 to be a local operator (e.g. at (x, —t ) in Eq. dlBli ). but this is 
not necessary for our discussion. Any interpolating operator for the pion on the chosen time 
slice would do equally well. 
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FIG. 5: Connected diagrams contributing at 0{a) contribution to the amplitude for the decay 

7T + —> l + Ul. 

Having determined Ao and hence the amplitude u Ue a (p U(: ) (M 0 ) a p v £ p(p £ ), the 0(a°) con¬ 
tribution to the decay width is readily obtained 

i-itree/ + . ,+ ^ &F \VJ? ft 2 A Y 

r o W ->■ Vvt) =--- m n m £ 1-~ . (20) 

87 t \ m~ J 

In this equation we use the label tree to denote the absence of electromagnetic effects since 
the subscript 0 here indicates that there are no photons in the final state. 


B. Calculation at 0{a) 

We now consider the one-photon exchange contributions to the decay n + l + v £ and 
show the corresponding six connected diagrams in Fig. [5] and the disconnected diagrams in 
Fig. [HI By “disconnected” here we mean that there is a sea-quark loop connected, as usual, 
to the remainder of the diagram by a photon and/or gluons (the presence of the gluons is 
implicit in the diagrams). The photon propagator in these diagrams in the Feynman gauge 
and in infinite (Euclidean) volume is given by 

e ik-(xi—X2) 

( 2^)4 p2 ■ ( 21 ) 

In a finite volume the momentum integration is replaced by a summation over the mo¬ 
menta which are allowed by the boundary conditions. For periodic boundary conditions, 
we can neglect the contributions from the zero-mode k — 0 since a very soft photon does 
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FIG. 6: Disconnected diagrams contributing at 0(a ) contribution to the amplitude for the decay 
7r + —>• £ + vi. The curly line represents the photon and a sum over quark flavours q, q\ and q 2 is to 
be performed. 


not resolve the structure of the pion and its effects cancel in T 0 — Tjf in Eq. (J3J). Although 
we evaluate r 0 + r 1 (A£ l ) (see Eq. (J2jl) in perturbation theory directly in inflnite volume, 
we note that the same cancellation would happen if one were to compute Ti(AE) also in a 
finite volume. Moreover from a spectral analysis we conclude that such a cancellation also 
occurs in the Euclidean correlators from which the different contributions to the decay rates 
are extracted. For this reason in the following r 0 and rjf are evaluated separately but using 
the following expression for the photon propagator in finite volume: 

2 _^ e ik-(xi-x 2 ) 

S^A(x ll X 2 ) = 5^ — 22 4 V cin 2 b. ’ 

k=^-n-,k^ 0 2 


where all quantities are in lattice units and the expression corresponds to the simplest lattice 
discretisation, k , n, x\ and x 2 are four component vectors and for illustration we have taken 
the temporal and spatial extents of the lattice to be the same (L). 

For other quantities, the presence of zero momentum excitations of the photon field is a 
subtle issue that has to be handled with some care. In the case of the hadron spectrum the 
problem has been studied in [22]] and, more recently in 3,[4|, where it has been shown, at 
0(a), that the quenching of zero momentum modes corresponds in the infinite-volume limit 
to the removal of sets of measure zero from the functional integral and that finite volume 
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effects are different for the different prescriptions. 

We now divide the discussion of the diagrams in Fig. [5] and Fig. [6] into three classes: those 
in which the photon is attached at both ends to the quarks (diagrams Eta)-|5](c) andHD(a), (b) 
and (d)), those in which the photon propagates between one of the quarks and the outgoing 
lepton (diagrams G3(e), [5](f ) and 13(c)) and finally diagram [5](d) which corresponds to the 
mass and wave-function normalisation of the charged lepton. We have already discussed the 
treatment of the wave function renormalisation of the lepton in detail in Sec. lIIIl so we now 
turn to the remaining diagrams. 

1. The evaluation of diagrams Fig. [ 5 |( a)-(c) and Fig.\B(a),(h) and (d) 

We start by considering the connected diagrams [5](a)-(c). For these diagrams, the leptonic 
contribution to the amplitude is contained in the factor [u Ue (p Ue ) 7"(1 — 7 5 )iy(p£)] and we 
need to compute the Euclidean hadronic correlation function 

C'i(t) = ~ I d 3 xd' L x l d 4 x 2 (0\T{Jw{0)jv(xi)jv,(x2)(f) ! (x,-t)}\0)A(x 1 ,X2). (23) 

where T represents time-ordering, Jfy is the V-A current dj u (l — y 5 ) u and we take —t < 0. 
jp is the hadronic component of the electromagnetic current and we find it convenient to 
include the charges of the quarks Qf in the definition of j: 

3u(x) = ^Qf f( x ) Infix), (24) 

/ 

where the sum is over all quark flavours /. The factor of 1/2 is the standard combinatorial 
one. 

The computations are performed in Euclidean space and in a finite-volume with the 
photon propagator A given in Eq. (1221) (or the corresponding expression for other lattice dis¬ 
cretisations). The absence of the zero mode in the photon propagator implies a gap between 
m n and the energies of the other eigenstates. Provided one can separate the contributions 
of these heavier states from that of the pion, one can perform the continuation of the cor¬ 
relation function in Eq. (I25|) from Minkowski to Euclidean space without encountering any 
singularities. From the correlation function C\ (t) we obtain the electromagnetic shift in the 
mass of the pion and also a contribution to the physical decay amplitude, as we now explain. 
For sufficiently large t the correlation function is dominated by the ground state, i.e. the 
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pion, and we have 

p—IUnt 

Co(t ) + C^t) ~ —- ^ (0 14,(0) 17T+), (25) 

where the electromagnetic terms are included in all factors (up to 0(a)). Writing = 
+ Sm n , where 5m n is the 0 (a) mass shift, 

e-mirt (26) 

so that C\(t) is of the schematic form 

Oi(t) = Co(t) (ci t + C 2 ). (27) 


By determining Ci we obtain the electromagnetic mass shift, = — ci, and from c 2 we 
obtain the electromagnetic correction to 20 (0 | Jty(O) 1 7 r + )/ 2 m 7r . Note that is gauge 
invariant and infrared finite, whereas the coefficient C 2 obtained from these diagrams is 
neither. 

In order to obtain the contribution to the 7 r —* ^ decay amplitude *4. we need to remove 
the factor (e~ mnt /2m n )Z^ on the right-hand side of Eq. (l25]h including the O(a) corrections 
to this factor. Having determined Ci, we are in a position to subtract the corrections present 
in rriif The 0(a) corrections to Z^ are determined in the standard way, by performing the 
corresponding calculation to C\(t) but with the axial current A replaced by </>: 

Cf(t) = ~ f d 3 xd A x 1 d A x 2 (0\T{(p(0,0)j fl (x 1 )j^x 2 )(j) ] (x,t)} | 0 ) A(xi,x 2 ) (28) 
= C^(t)( Cl t + cf). (29) 


Cl 


We finally obtain 

z* = zs(i + t^-^ 

and the 0 (a) contribution to the amplitude from these three diagrams is 


(30) 


5A 


Ao 



—)■ 

2 ) 


(31) 


For these three diagrams the 0(a) term can be simply considered as a correction to f n . 
Note however, that such an would not be a physical quantity as it contains infrared 
divergences. 

The treatment of the disconnected diagrams in Figs.| 6 ](a), (b) and (d) follows in exactly 
the same way. These diagrams contribute to the electromagnetic corrections to both the pion 
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mass and the decay amplitude in an analogous way to the discussion of the connected dia¬ 
grams above . It is standard and straightforward to write down the corresponding correlation 
functions in terms of quark propagators. We do not discuss here the different possibilities 
for generating the necessary quark propagators to evaluate the diagrams; for example we 
can imagine using sequential propagators or some techniques to generate all-to-all quark 
propagators. 


2. The evaluation of diagrams Fig.\$(e)-(f) 

For these diagrams the leptonic and hadronic contributions do not factorise and indeed the 
contribution cannot be written simply in terms of the parameter f n . We start by considering 
the Minkowski space quantity 

u Vta (p Vl )(Mi) a pVei 3 (pi) = - jd A xi d A x 2 (0| T^aq) J^(0)) |7r) (32) 

x iD M (x 1 ,x 2 ){u Vt (p l , l )'f (1 - l b ){iS M {x 2 ))Yvi{Pi)}e m ' X2 , 

where iSm and iD M are the lepton and (Feynman gauge) photon propagators respectively 
in Minkowski space (more precisely the photon propagator with Lorentz indices (p, a) is 
iD m g p a i but the Lorentz indices have been contracted with the electromagnetic currents 
in ((32])). In order to demonstrate that we can obtain the 0(a) corrections to the decay 
amplitude from a Euclidean space correlation function, we use the reduction formula to 
rewrite the expression in Eq. (132|) as 

u^aipvJiMjapVepipi) = i lim (k 0 2 - w 2 ) f d 4 x 1 d i x 2 d i x e~ lk ° x ° 

ko ym n J 

(0|T(^(xi) J^(0)7r(x))|0) iD M (xi,x 2 ) [u Vl (p Vl )^ v (l - r y 5 )(iS M (x 2 ))Yve(pe)\e m ' X2 , (33) 

where n(x) is the field which creates a pion with amplitude 1. On the other hand the 
Euclidean space correlation function which we propose to compute is 

Ci(t) a p = - l<l 3 xd i x 1 d i x 2 <0|T{J^,(0) -f)}| 0) A(x 2 ,x 2 ) 

Here S and A are Euclidean propagators, and a,ff are spinor indices. Similarly to the 
discussion in Sec. IIVB 11 provided that the pion is the lightest hadronic state then for large 
t, Ci(t) is dominated by the matrix element with a single pion in the initial state. 
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FIG. 7: Zoom of the lepton-photon vertex at x 2 from the diagrams in Fig.[5|e) and (f). 


In view of the factor e Ett2 on the right-hand side of Eq. (1 341) . the new feature in the 
evaluation of the diagrams in Fig.[5](e) and (f) is that we need to ensure that the £2 integration 
converges as ^ 2 1 —>■ 00 . For £2 < 0 the convergence of the integral is improved by the presence 
of the exponential factor and so we limit the discussion to the case t 2 — !> 00 . Eg = s/rrij + p e 2 
is the energy of the outgoing charged lepton with three-momentum pg. To determine the 
t 2 — > 00 behaviour, consider the lepton-photon vertex at x 2 from the diagrams in Fig.^e) 
and (f), redrawn in Fig.[7| kg and k y are the four-momentum variables in the Fourier 
transform of the propagators S(x 2 ) and A(xi,x 2 ) respectively in Eqs. (1321) - (1341) . The t 2 
integration is indeed convergent as we now show explicitly. 

1. The integration over x 2 implies three-momentum conservation at this vertex so that in 
the sum over the momenta kg + k^ — pi, where pg is the momentum of the outgoing charged 
lepton. 

2. The integrations over the energies k^g and /c 4 7 lead to the exponential factor e“( u ’ i+a '' 1, ) t2 , 
where uig = yjk t 2 + m 2 , u; 7 = \jk^ + m*, and m 7 is the mass of the photon introduced as 
an infra-red cut-off. The large t 2 behaviour is therefore given by the factor e -(<^+ w 7 - E d *2 _ 

3. A simple kinematical exercise shows that in the sum over k 1 (with kg = pi — fc 7 ), the 
minimum value of u)g + tn 7 is given by 


(ui + w 7 ) ml „ = \J (m, + m 7 ) 2 + Vi ■ 


(35) 


4. Thus for non-zero m 7 , the exponent in f or l ar g e t 2 j s negative for every term 

in the summation over k 1 and the integral over t 2 is convergent so that the continuation 
from Minkowski to Euclidean space can be performed. 

5. We note that the integration over t 2 is also convergent if we set m 7 = 0 but remove the 
k = 0 mode in finite volume. In this case cog + cn 7 > Ei + [1 — (pe/Eg)]\k m i n \. 

In summary the t 2 integration is convergent because for every term in the sum over 
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momenta ay + ay, > Ei and so for sufficiently large t we can write 

( 36 ) 

and the contribution from the diagrams of Fig.EJ^e) and Eff) is u a (p Ve )(Mi) a pVp(p£). This 
completes the demonstration that the Minkowski-space amplitude (133|) is equal to the pion 
contribution to the Euclidean correlation function (l34jh up to a factor Zq which accounts 
for the normalisation of the pion held. 

Again the evaluation of the correction to the amplitude from the disconnected diagram 
in Fig. Etc) follows in an analogous way. 


V. CALCULATION OF r pt (A£’) 


The evaluation in perturbation theory of the total width T pt = in infinite volume, 

was performed by Berman and Kinoshita in 1958/9 If], 2^|, using the Pauli-Villars regulator 
for the ultraviolet divergences and a photon mass to regulate the infrared divergences in both 
rf and rf. rf is the rate for process 7r + —y for a pointlike pion with the energy of 

the photon integrated over the full kinematic range. We have added the label pt in r pt to 
remind us that the integration includes contributions from regions of phase space in which 
the photon is not sufficiently soft for the structure of the pion to be reliably neglected. We 
do not include this label when writing Ti(A.E) because we envisage that A E is sufficiently 
small so that the pointlike approximation reproduces the full calculation. 

In our calculation, Tp is evaluated in the W-regularisation, so that the ultra-violet diver¬ 
gences are replaced by logarithms of M w . For convenience we rewrite here the expression 
for r pt (A£) from Eq. (fTfil) 


r pt (A£) = r pt + t ^ ae ) = r* ree + r“’ pt + r[, d) ’ pt + t ^ ae ) . ( 37 ) 

r^ 66 and r{, d) ’ pt have already been presented in Eqs. (1201) and (1T5|l respectively. In the 
following we give separately the results of the remaining contributions to r pt (Ai7) also 
using a photon mass m 7 as the infrared regulator. We neglect powers of m 7 in all the 
results. 

In the perturbative calculation we use the following Lagrangian for the interaction of a 
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and 



FIG. 8: One loop diagrams contributing to the wave-function renormalisation of a point-like pion. 


point-like pion with the leptons: 


= i G F firV* d {(df, - ieA^n} <ip, 


1 + 75 


Vl 2 


'Y^'ipe > + Hermitian conjugate. (38) 


The corresponding Feynman rules are: 


7r 1 



t + 


vi 


-iG F f„v: d ri 


(39) 


7r 


+ 



l + 


Vl 


ieG F f,V: d g“" 


In addition we have used the standard Feynman rules of scalar electromagnetism for the 
interactions of charged pions in an electromagnetic field. 

We start by giving the 0(a ) contributions to Fo’ pt . 

• Wave function renormalisation of the pion: The contribution of the pion wave function 
renormalisation to Fq' p1; is obtained from the diagrams in Fig. [8] and is given by 

U = r where Z, = -2 log (|J) - 2 log (g) - |. (40) 

These diagrams correspond to those in Fig. [5](a), Fig. [5](b) and Fig. [5](c) in the composite 
case. 

• 7 t- i Vertex: The remaining graphs contributing to Fg’ pt are the n - £ vertex corrections 
from the diagrams shown in Fig. [9] and their complex conjugates. The contribution from 
these diagrams is 


a 


r[p = r“ ee x — Zv-i where 


47T 


(41) 


/ 7 

/ mz 


mz 


z, = - 2 I±7 log (rf) log ( ^ J + 4 log v ^ 


1 - ri 


+ 


log 2 (r, 2 ) + 2++h log (rj) - 1 

j -Id -L I D 


( 42 ) 
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FIG. 9: Radiative corrections to the pion-lepton vertex. The diagrams represent 0(a) contribu¬ 
tions to Tq 1 . The left part of each diagram represents a contribution to the amplitude and the 
right part the tree-level contribution to the hermitian conjugate of the amplitude. The correspond¬ 
ing diagrams containing the radiative correction on the right-hand side of each diagram are also 
included. 








FIG. 10: Diagrams contributing to Ti(A.E). For diagrams (c), (d) and (e) the “conjugate” con¬ 
tributions in which the photon vertices on the left and right of each diagram are interchanged are 
also to be included. 

and rg = These diagrams correspond to the diagrams Fig. [51(e) and Fig. [51(f) in the 

composite pion case. 

Next we give the contributions to Ti(A.E) where the real photon is emitted and absorbed 
by the pion (mr), the charged lepton (it) or emitted by the pion and absorbed by the lepton 
or vice-versa (wt). The results are presented in the Feynman gauge: 

^2^(k,r)£ v (k,r) = g„ v , (43) 

r 

where e M (fc,r) are the polarisation vectors of the real photon carrying a momentum k, with 
k 2 = 0 in Minkowski space. 

• Real photon emission, m r: The contribution to r 1 (AFl) from the emission and absorption 
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of a real photon from the pion, represented by diagram (a) in Fig.^Ol is given by 

ry 

= r|, ree x — (RY + RY) , where 


(44) 


RY = 2 log 


m 2 \ 
4A E 2 ) 


+ 4, 


r>7T7r 

/l 2 


2^ 4 


(1 


-2W 


log(l - r E ) + 


f e (6 -r E - 4r 2 ) 


(1 - r?)* 


(45) 


v E = 2A E/m n and 0 < re < 1 — r 2 . Here we have separated RY, the contribution in the 
eikonal approximation from RY which vanishes as A E —y 0. In the eikonal approximation 
only the leading terms in the photon’s momenta are kept in the numerator and denominator 
of the integrand as r E —> 0. RY contains the infrared divergence. 

• Real photon emission, it. The contribution to Ti(AE) from the emission and absorption 
of a real photon from the charged lepton, represented by the diagram (b) in Fig. flOl is given 
by 


rf = YY e X (R.f + R-f) , where 


(46) 


R 


tt 

i 


2 log 


m: 


_A) 

4A E 2 ) 


- 2 


1 + r. 


log (A 2 ) 


and 


R 


it 

2 


1 + (4 r E - 6)r 2 
(1 - r|) 2 


log(l - r E ) - 


r E ( 6 — 3r£ — 20r|) 
f 2(1 - r 2 ) 2 


r E (r E + 4 rf) 
(1 - r?)* 


log(r^) 


(47) 


• Real photon emission, ni: Finally, the contribution to Ti(AE) from the emission of a real 
photon from the pion and its absorption by the charged lepton, represented by the diagrams 
(c) - (f) in Fig. [TOJ is given by 


r 


7 rl 
1 


rr x x «+ Rf ), 


(48) 


where 


< = 2i±I|log(r ( 2 )log 


' 1 - V t 


m. 


\ 1 + rj 


AAE 2 




l ~ r t 


[log(' 


2,1 2 ,1 + ri T . ^ 2 , . 

- ^ -4--1 Li 2 (1 — ri) and 

1 - r i 


tdtt£ _ 

— 


2 r E + rj — 2 


(1 


-2N2 


log(l - 7 \e) + 


4r ; 


(1 


~ 2\2 


log(r £ 2 ) + 


r E ( 2 + r E ) 
(1 - r?)=* 


1 + r? 


Li 2 (^e) • 
(49) 
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Note that for diagrams (c), (d) and (e) we include the conjugate contribution in which the 
photon vertices are interchanged between the left and right parts of the diagrams. Thus for 
example, in addition to diagram (c) there is the diagram in which the photon is emitted 
from the lepton on the left and absorbed on the pion on the right. 

We are now in a position to combine the results in Eqs. (1401) - (]49l) in order to obtain the 
final expression for T pt (Aif). As expected the infrared cutoff cancels and we find 

’ 2 x 2 - 10r 2 


r pt (A£) = rr x ( 1 + <j Slog ) + log (rf) - 41og(r|) + 


1 - r i 


■l°g(r £ ) 


„ 1 + r 2 n 2 2 ,1 + rj 

-2--4 log(r|) log(r 2 ) - 4- 


1 — ri 


1 — ri 


Li2(l 


3 + r| - 6rf + 4r E (-l + rf) r E (4 - r E 

- log(l -r E ) + 


4 77 


(1 - r 2 ) 2 


(1 - r 2 ) 2 


log {rl 


r E (—22 + 3r E + 28r 2 ) 1 + r 2 . , . 

— 4--4 Li 2 (r E ) 


T-r 2 


(50) 


2(1 - r 2 ) 2 

Note that the terms in square brackets in eq. (150]) vanish when r E goes to zero; in this limit 
T pt (A£) is given by its eikonal approximation. 

The total rate is readily computed by setting r E to its maximum value, namely r E = 1—r 2 , 
giving 

' 2 x 3 rj 


r pt = r 


tree 

0 


x<1+ s ( 31og (w 


- 81og(l - rf) - 


(1 - rjy 


log(r £ ) 


(51) 


■'^4 Li 2 (l 


r]) + 


13 — 19r 2 6 — 14r 2 — 4(1 + r 2 ) log(l — r 2 ) 


+ 


log(r 



±-'e 2(1-r 2 ) 1-r 2 

The result in Eq. (T5TT) agrees with the well known results in literature jl6, 17], which provides 
an important check of our calculation. We believe that the result in Eq. (1501) is new. 

In the description of our method above, we limit the photon’s energy to be smaller than 
A E to ensure that the photon is sufficiently soft for the pointlike approximation to be valid 
in the evaluation of Ti(A12). It is of course possible instead to impose a cut-off on the energy 
of the final-state lepton, requiring it to be close to its maximum value Ef dX = ^(1 + r 2 ). 
For completeness we also give, up to O(AEe), the distribution for T pt (AlA) defined as 

dT pt 


_ j^max 


T pt (AE e ) = I dE' 

J E ™ ax -AE e 


dE'i ’ 


( 52 ) 
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where 0 < A Eg < (m n — mf) 2 /(2m ff ); 

r-*(AB ( ) = +” x 1 1 + £■ 3 log + 8 log (1 - rf) - 7 

3-7 'rf + 8A£> + 4(1 + r^) log (1 - r^) 


log {rj) 


1 - r I 


log (2AEg) ( -8 - 4 1 + ^ log (rf) 


l~ r l 


(53) 


VI. REGULARISATION AND CANCELLATION OF INFRARED 
DIVERGENCES IN FINITE-VOLUMES 


In the previous section we have explicitly demonstrated the cancellation of infrared di¬ 
vergences in the perturbative quantity r pt (Alf). This of course is simply the standard 


Bloch-Nordsieck cancellation 
of infrared divergences in 


ll| . In this section we discuss in more detail the cancellation 


at 0 (l) = f q — r 


-iQ,pt 


(54) 


We have already explained in Sec. lIIIl that the contribution of the lepton’s wave function 
renormalisation in ATo(L) is simply to introduce the tilde in Tq, denoting that the corre¬ 
sponding contribution to the matching factor between the lattice and W -regularisations is 
to be removed. We also do not discuss further the evaluation of the remaining infrared-finite 
terms in the matching factor because these are straightforward to evaluate (see e.g. Eq. (fT4j) 
for the Wilson action). Here we concentrate on the remaining diagrams in Figs. [5] and [6] and 
the corresponding diagrams for the point-like meson. 

Although the right-hand side of Eq. (l54li is a difference of decay widths, since at this 
order the widths are linear in the 0(a) virtual amplitude, we can equivalently consider 
the difference of the 0(a) contributions to the amplitudes. In order to reduce statistical 
fluctuations when performing the sum over the gauge field configurations, we define the 
ratios 


R a 


A a 



R a,pt 


yp.pt 



(55) 


where A a and A a,pt are the 0(a) amplitudes corresponding to the widths in Eq. (154|) . The 
non-perturbative amplitude A a is precisely the quantity that we propose to compute nu- 
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merically in a lattice simulation. It is then combined with A a,pt , for which we have given 
the explicit expression in infinite volume in Sec.lVl 

In the calculation of A a,pt we set the mass of the photon to zero and consider the theory 
on a finite volume of length L, which will be used as an infrared regulator. The form of 
the vertices and propagators is the same as in the infinite volume (the ultraviolet cutoff 
is provided by the W-regularisation), but the momenta are quantized k M = 2ir/L x = 
27r/ (Na) x n M where — oo < < +oo and N is the number of lattice sites in one direction, 

which for simplicity we take to be the same in all directions. 

The calculation of Aq is performed non perturbatively on the same finite volume as in the 
perturbative case, but in a numerical simulation and with the photon propagator defined as 
in Eq. (|22l) . which does not contain the zero mode. Indeed, as already discussed in Sec. IIV B1 
the zero mode does not contribute to the difference 

A R(L) = R a - R a ’ pt . (56) 

This is a gauge invariant, ultraviolet and infrared finite quantity and for these reasons we 
expect that its finite volume effects are comparable to those affecting the 0(a) corrections 
to the hadron masses (that are also gauge invariant, ultraviolet and infrared finite). The 
formalism introduced in this paper was necessary because To and Ti are separately infrared 
divergent. 

We should add that in principle any consistent regularisation of the infrared divergences 
is acceptable. The main criterion for the choice of the infrared regulator will be determined 
by the precision of the terms remaining after the cancellation of the infrared divergences in 
a numerical simulation. 


VII. SUMMARY AND PROSPECTS 

Lattice calculations of some hadronic quantities are already approaching (or even reach¬ 
ing) 0(1%) precision and we can confidently expect that the uncertainties will continue to be 
reduced in future simulations. At this level of precision, isospin-breaking effects, including 
electromagnetic corrections, must be included in the determination of the relevant physical 
quantities. In this paper we present, for the first time, a method to compute electromag¬ 
netic effects in hadronic processes. For these quantities the presence of infrared divergences 
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in the intermediate stages of the calculation makes the procedure much more complicated 
than is the case for the hadronic spectrum, for which calculations in several different ap¬ 
proaches (3-8] already exist. I11 order to obtain physical decay widths (or cross sections) 
diagrams containing virtual photons must be combined with those corresponding to the 
emission of real photons. Only in this way are the infrared divergences cancelled. We stress 
that it is not sufficient simply to add the electromagnetic interaction to the quark action 
because, for any given process, the contributions corresponding to different numbers of real 
photons must be evaluated separately. 

We have discussed in detail a specific case, namely the 0(a) radiative corrections to the 
leptonic decay of charged pseudoscalar mesons. The method can however, be extended to 
many other processes, for example to semileptonic decays. The condition for the applicability 
of our strategy is that there is a mass gap between the decaying particle and the intermediate 
states generated by the emission of the photon, so that all of these states have higher energies 
than the mass of the initial hadron (in the rest frame of the initial hadron). 

In the present paper, we have limited the discussion to real photons with energies which 
are much smaller than the QCD scale Aqcd- This is not a limitation of our method and 
in the future one can envisage numerical simulations of contributions to the inclusive width 
from the emission of real photons with energies which do resolve the structure of the initial 
hadron. Such calculations can be performed in Euclidean space under the same conditions 
as above, i.e. providing that there is a mass gap. 

In the calculation of electromagnetic corrections a general issue concerns finite-size ef¬ 
fects. In this respect, our method reduces to the calculation of infrared-finite, gauge-invariant 
quantities for which we expect the finite-size corrections to be comparable to those encoun¬ 
tered in the computation of the spectrum. This expectation will be checked in forthcoming 
numerical studies and studied theoretically in chiral perturbation theory. Indeed an analyt¬ 
ical calculation of the Unite-volume effects requires a detailed analysis of the form factors 
parametrising the structure dependent contributions (see. Eq. ( 1 B 4 D ). 

Although the implementation of our method is challenging, it is within reach of present 
lattice technology particularly as the relative precision necessary to make the results phe¬ 
nomenologically interesting is not exceedingly high. Since the effects we are calculating are, 
in general, of 0(1%), calculating the electromagnetic corrections to a precision of 20% or 
so would already be more than sufficient. As the techniques improve and computational 
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resources increase, the determination of both the QCD and QED effects will become even 
more precise. We now look forward to implementing the method described in this paper in 
an actual numerical simulation. 
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Appendix A: Matching between Lattice and W-regularisation 

In this appendix we briefly describe the matching between the lattice and W regularisa- 
tions, in perturbative QED for the complete basis of four-fermion operators 

Oxy = (dT xU ) (P iT y t) = T x ® T y , (Al) 

where Tx,y are Dirac matrices. We consider the following basis of five four-fermion operators 
given in Eq. (ITT]) : 

Oi = 7 m (1 - 7 5 ) (g) 7 m (1 - 7 5 ), 0 2 = + 7 5 ) < 8 > 7 M (1 - 7 5 ), 

0 3 = (1 - 7 5 ) <g) (1 + 7 5 ), 0 4 = (1 + 7 5 )® (1 + 7 5 ), ( A2 ) 

0 5 = ( 7^(1 + 7 5 ) < 8 ) 0-^(1 + 7 5 ). 

The complete basis is made up of ten operators. The five additional operators are obtained 
from ()\ -O 5 by the exchange (1 — 7 s ) tA (1 + 7 s ). Since the neutrino is electrically neutral 
its chirality is conserved and the operators Oi - O 5 do not mix under renormalisation with 
the remaining 5 operators and invariance under parity transformations ensures that the two 
5x5 renormalisation matrices are equal. For this reason, in the following we focus the 
discussion on the five operators of Eq. (1A2I) . Moreover, the basis of operators in Eq. lA2l is 
the complete basis of operators for a left-handed neutrino. 

With regularisations which respect chiral symmetry the four-fermion operator relevant 
for the leptonic weak decay, Oi, renormalizes multiplicatively. I 11 this appendix we are using 
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the lattice theory with Wilson fermions to illustrate the matching between the lattice and 
W-regularisations and the explicit breaking of chiral symmetry with this discretisation of 
QCD leads to the mixing of O\ with the other four operators O 2 -O 5 . If instead of using 
Wilson fermions, we used a lattice formulation with good chiral properties, such as domain 
wall fermions, the corresponding discussion to the one presented below would be restricted 
to the single operator 0\ which transforms as the (8,1) representation under SU(3 )lxSU(3)r 
chiral symmetry for the quarks. 

We define Zij(aM\y) to be the matrix which relates the operators 0* (i=l-5) in the 
lattice and W-regularisations: 


Of{M w ) = Zij(aM w ) 0^\a ). 


(A3) 


In order to perform the matching we adapt the RI-MOM renormalisation procedure devel¬ 
oped for QCD [24], although, as described below, all the calculations here are performed in 


perturbation theory. Let A) att and A) 1 (z=l-5) be the amputated 4-quark Green function 
of the operator O t with the lattice and W regularisations respectively, both with external 
momenta p as illustrated in Fig. [IT] We determine Z by imposing that 


(zd zp zp) Z lk Tr (A f“Pj) = Tr (A f P,) . (A4) 

The projectors Pj are dehned by their action on the tree-level Green function A-° , 

Tr (AfP,) = <% , (A5) 

where the trace here and in Eq. (1A4D is dehned by Tr(A,Pj) = Tr {Y l x P x Fy Py) for 

Oi = Tj^ <g> Ty and Pj = P 3 X <S> Py. Z u ^ are the matching factors for the wave function 

1 . 

renormalisation constants of the corresponding fermion fields, e.g. = Z£ u latt . 

Consider the perturbative expansion of the amputated bare Green function in powers of 
the electromagnetic coupling in either the lattice or W-regularisations, 


A,=A ‘ 0) + s A<w+ -- (A6) 

In order to implement the matching conditions between the two regularisation schemes we 
require the quantities Tr (AjPQ in both schemes. At one-loop order we write 

Tr (A‘ att( 1 ) P ? ) = Dij and Tr (kf (V) P^j = C tJ . 


(AT) 
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FIG. 11: One-loop Feynman diagrams computed for the renornralisation of the four-fermion oper¬ 
ators Oxy = (dTx u ) (iy Ty t) = Tx <8> Fy. 


We represent the matching of the wave functions in the lattice and W -regularisations up 
to one-loop order by Z q = 1 + (a/47r) Zp + .... Using Eq. (1A4|) . we see that the matching 
matrix of the operators in Eq. (1A2h at 0(a) is given by 

4 1 ’ = C ‘i - D ‘i + \ (4 1 * + Z™ + zf) <% . (AS) 

We have presented the O(a) contribution to the matching factor for the wave function 
of the charged lepton in Eq. (1T5D of Sec. lIIIl 


Zp ] = -3/2 - log (a 2 Mp) - 11.852 , 


(A9) 


zP and ZP differ from Zp only by factors of Q ^ and Q 2 d , where Qf is the charge of 
the fermion /. We have verified with an explicit calculation that the contribution to the 
matching given in eq.(A9) is the same whether evaluated for an on-shell or an off-shell 
external lepton. 

In order to evaluate the matrices C i3 and Dij it is necessary to compute the Feynman 
diagrams shown in Fig. [II] in the two regularisation schemes. All the external momenta are 
chosen to be equal to p and all external particles are taken to be massless. We deduce 
from the results of the corresponding QCD calculation performed in j^. (Ref. 25] includes 
a package containing an ASCII hie, in order to make the results most easily accessible to 
the reader.) Diagrams 1, 2, 3 of Fig. [Til correspond to the diagrams d 5 , d e and d\ of 25 ]. 


The expression for the lattice wave function renormalisation can be obtained from 


26]. 


We now present results for the standard Wilson fermions and the “naive” QED gauge 
action, for which the tree-level lattice photon propagator in the Feynman gauge is given 
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in Eq. (12 2 p . In infinite volume the sum over momenta in Eq. (122|) is replaced by the corre¬ 
sponding integral. By combining the ingredients discussed above, we obtain the following 
result for the 0(a) contribution to the renormalisation matrix of Eq. (1A8I) : 

/ 2 L 


Z = 


W - 15.539 

0.536 

1.607 

-3.214 

-0.804 

0.536 

L w - 14.850 

-3.214 

1.607 

-0.402 

0.402 

-0.804 

-§ L w - 13.702 

-1.071 

0 

-0.804 

0.402 

-1.071 

-f L w - 13.702 

72 Lw — 0.( 

-9.643 

-4.822 

0 

AL W - 2.756 

f Lw ~ 15. 




(A10) 

where L w = log (a 2 M^). 

The four-fermion operator relevant for the leptonic decay rate is 0±. From Eq. (1A10D we 
obtain the expression in Eq. (fTUj) for 0\ in the VE-regularisation in terms of the bare lattice 
operators. 

The result presented for Z in Eq. (1A10I) above is also valid if the twisted-mass (or 

n 

Osterweilder-Seiler [271]) lattice regularisation is used for the fermions instead of the Wilson 
action. This statement follows from the observation that the twisted mass action, in the 
so called twisted basis j^, only differs from the Wilson action by the presence of y 5 in the 
mass term. The two actions are therefore identical in the chiral limit and all renormalisa- 
tion constants are equal for Wilson and twisted-mass fermions in the twisted basis in all 
mass-independent renormalisation schemes. The renormalisation constants for twisted-mass 


fermions in the p 
twisted rotation 


rysical basis are obtained from those in the twisted basis through a simple 


The lattice results in 


25 


26] are also given for a number of pure gauge actions including 


the tree-level Symanzik and Iwasaki actions. For completeness we give below the results for 
the renormalisation matrix for these two choices of the gauge action: 

\ 


7(1) _ 

Z / T5 — 


'\y — 12.399 

0.451 

1.354 

-2.709 

-0.677 

0.451 

L w - 11.866 

-2.709 

1.354 

-0.339 

0.339 

-0.677 

-\L W - 10.978 

-0.903 

0 

-0.677 

0.339 

-0.903 

-§ L w - 10.978 

j^Lw — 0.( 

-8.127 

-4.063 

0 

\L W - 2.132 

fAv-12. 


(All) 
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7(1) _ 

A Iw — 


2 L 


W - 11-732 

0.323 

0.969 

-1.938 

-0.485 

0.323 

L w - 11.525 

-1.938 

0.969 

-0.242 

0.242 

-0.485 

— — 11.181 

-0.646 

0 

-0.485 

0.242 

-0.646 

— — 11.181 

toR 

R 

1 

0 

-5.815 

-2.908 

0 

4 L\y — 0.826 

20 T 11 

-9 ~ 


\ 


(A12) 


Appendix B: Structure dependent contributions to the real decay 


In this appendix we estimate the size of the neglected structure-dependent contributions 


to the decay P + —> l+Vf.'y for light mesons, P + = ir + ,K + . We 


3 ase our estimates on 


the results of the phenomenological analyses performed in Refs. [29M31I] based on the use 
of chiral perturbation theory at 0(p 4 ). Although the relevant expressions have also been 
derived at 0(p 6 ) [32,133] (see also page 10 of |34j), in that case there are too many unknown 


low-energy constants to be useful in making an estimate. As was done in the main body 
of the paper, for the general framework we give the explicit formulae for pion decays; the 
generalisation of the framework to kaons, and indeed also to H-mesons and R-mesons decays 
is straightforward. We then make the numerical estimates of the structure dependent effects 
for pions and kaons based on chiral perturbation theory. Finally we make some comments 
about structure dependent terms when P + is a heavy-light meson, D + or B + . 

The starting point of the analysis is the decomposition in terms of Lorenz invariant form 
factors of the hadronic matrix element (see also Eq. 


H^(k, Pn ) = I d A xe ikx T(0\f(x)J^(0)\ir( Pn )) . 


(Bl) 


We follow the standard convention of separating the contribution corresponding to the ap¬ 
proximation of a point-like pion (also frequently called inner bremsstrahlung) , from the 

structure dependent part H 

(B2) 


H^ = H& + Hg 


Hp t is simply given by 


H^ - f 

pt Jtt 


9^- 


(2 Pn - kY( Pn - ky 
(Ptt - k ) 2 - ml 


(B3) 
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The structure dependent component can be parametrised by four independent invariant 
form factors which we define as 


#SD = H, [k 2 g^ - k»k v ] + H 2 { [{k • p n - k 2 )W - k 2 (p n - ky] (; V, - k) u } 

Fa 


.F v 


-i— e^k aPnp + — [(A; ■ p. - k 2 )g^ - ( Pn - kyk u ] 

uItt h m n 


(B4) 


Note that the vector Ward Identity k M H^ v = f n p", derived in Ref. [29|, is saturated by 

k,H£ = y P :, k,HZ = 0. (B5) 


As discussed in the main body of the paper, also contains the infrared divergences which 
appear in the virtual- and real-photon contributions to the decay rate. These observations 
motivate the decomposition in Eq. (11121) . 

In the calculation of the decay rate for n + —> £ + uff the tensor is contracted with 
the polarisation vector of the real photon. In physical gauges with e* ■ k = 0 we define 


H v = e*H^ , 


(B6) 


so that 


H v S d = 


i— e^k a p^ - — [(k ■ pi, - k 2 )g^ - (p v - k)^] \ , (B7) 


m, ^ m 

showing that the structure dependent part of the decay rate can be parametrized in terms 
of the two form factors Fy and Fa- 

Before performing the integrations over the three-body phase space, the differential decay 
rate can be expressed as a function of the two independent Dalitz variables ( p n = pi+p v + k) 

(Pn ~ k) 2 


_ ( Pn-Pe ) 2 1 

X( ~ - o - r i , 


x 7 = 


+ 1 


(B8) 


m* m‘ 

The decay rate as a function of the photon’s energy in the pion’s rest frame can be obtained 
by performing the integration over xi with the limits Xi G [x™ m , x™ ax ] where 

1 — x^ — r 2 p 


xT n = 


2(1-x 7 ) 


x 1 -r 2j r J (x 7 + r 2 ) 2 - 4r 2 


1 _ rv* _ ip" 

max _ i 2 1 x 7 ' i 

x e ~ 1 '7 


x 1 -r - J (x 7 + r 2 ) 2 - 4 r 2 


(B9) 


2(1 — Xj) 

r£ = mi/niTr and r 7 = m 7 /m 7r . The total decay rate is obtained by performing the integral 
over x 7 in the range x 1 G [x™ m , with 

2 


OX 


= r 7 (2 -r 7 ) , 


x „ 


= 1 — To 


(BIO) 
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The photon’s mass m 1 was introduced in the definition of r 7 to regulate the infrared di¬ 
vergences in the point-like contribution. For the structure dependent contribution, which 
is infrared finite we can set m 7 —y 0 and simplify the above expressions by making the 
replacements 


xf m ^ (1 - x 7 ) + 


2 


(1 


ay 


Xo 


^ 1 


ay 


i —y 0 , 


ay 


i-)- 1 — r 


i ■ 


(BH) 


The different contributions to the differential decay rate have been obtained in Ref. 29]. 
Writing T [ = T^ + T^ D + T* nt , where T™ t is the contribution to the decay rate coming 
from the interference between the point-like and the structure-dependent amplitudes, we 
confirm the following results: 

4?r d 2 Tf _ 2f pb (xy,x t ) 
a Tg ree dx-fdxi (1 — rf) 2 

4tt d 2 Tf° _ ml {[F v (xy) + F A (xy)] 2 f£ D (xy,xe) + [F v (xy) - F A (x 1 )} 2 fg D (x 1: x e )} 


a To ree dx^dxi 

4tt d 2 Tf T 
a Tq 1 " 66 dx^dxa 


2f 2 r 2 (l-r 2 Y 

2m n { [Fv(xj) + F A (xy)\ f^ T (xj, ay) + [Fy(a: 7 ) - F A (xy)\ T (x 7 , ay)} 


U (1 - rj) 2 


(B12) 


The functions appearing in Eq. (1B12|) are 

1 - ay 


fpt(xj, xe) = 


X; 


l(Xry + X£~1) 


x 7 + 2(1 -z 7 )(1 -r e ) 


2 \ 2 x^(1-rl) 


Xry + X£ — 1 


/ s h D (x 7 , Xg) = (x 7 + X£-1) [(x 7 + Xe - 1 + ry 2 )( 1 - x 7 ) - rj\ , 

/sd( x 7 , x z) = “(! ~xe) [(ay - 1 + r 2 )(l - x 7 ) - r 2 \ , 

l- x e 


(B13) 


«/*INt(^7’ Xi) ~ 




X. 


,(x 7 + X£~ 1) 


\(xy + ay - 1 + r 2 )(l - Xy) - r 2 ] , 


xe 


,{xy + ay - 1) 


x. 


[x 2 + (xy + xe - 1 + rf)(l - Xy) - r 2 ] 


Whilst we confirm the results of Ref. 29J, we note t 


no 

e t 

35, 

37 


rat we disagree with the sign of the 


The sum of Eqs. (I44|h (1461) and fl48l) of the main body of the paper can also be obtained 
by integrating the point-like contributions over ay with the limits given in Eq. (1B9[) and over 
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x 1 in the range [r 7 ( 2 — r^),r E \. It will be useful below to define the following quantities, 


Ql( x -y) = 


Rf(AE) = 


47t dT^(x^ 
dx~/ ’ 

rf(AE) 


a T^ 66 


pa,pt 

1 0 


+ r { Q ] ’ pt + rf (A e) 

where A E = 777^/2 and Tg’ pt and TQ^’ pt have been dehned in the main body of the paper 
(see Eq. (13T1) ). Notice that the quantity in the denominator of Rf(AE) is infrared finite 
(although it does depend on M w , the ultraviolet cutoff in the VE-regularisation). 

In the following we use phenomenological parametrisations of the form factors in order to 
estimate the size of the structure-dependent contributions to the decay rate IY For the case 
of light mesons, we can use the results of the calculations of refs. 29-31] (see also ref. 0) 
based on chiral perturbation theory and approximate the form factors as constants. At 
0(p A ) in chiral perturbation theory, 

rrip 8m p 


A = {pt,SD,INT} , 
A = {SD,INT} , 


(B14) 

(B15) 


Fv = 


and Fa = 


(L r 9 + L 


lot > 


(B16) 


4vr 2 f n U 

where P = 7r or K and Lg, L[ 0 are Gasser-Leutwyler coefficients. The numerical values 
of these constants have been taken from the review by M.Bychkov and G.D’Ambrosio in 
Ref. 35]; the values of Fy and Fa are 0.0254 and 0.0119 for the pion and 0.096 and 0.042 for 


the Kaon (for the pion these values o 
can be found in the supplement to 


le form fact ors, obtained from direct measurements, 


351 ] found in 36]). In Figs.fl2l and fl3l we compare 


the point-like, structure-dependent and interference contributions to the decays 7r —> 
and K — y fz/y respectively. As can be seen, interference contributions are negligible in 
all the decays. The structure-dependent contributions can be sizeable because they are 
chirally enhanced with respect to the point-like contribution (notice the factor 1/rf in the 
second equation in (IB 121) ). From the phenomenological estimates of the form factors, this 
happens for the real decay K — > eryy. On the other hand, for < 20 MeV both structure 
dependent and interference contributions can be safely neglected with respect to the point¬ 
like contributions for all the decays of pions and the decay K —> /rzzy. We learn from 
Refs. jl4, 15] that a cutoff on the energy of the photon in the rest frame of the decaying 
particle of 0(20 MeV) is experimentally accessible. 

The application of chiral perturbation theory described above does not apply to the decays 
of D and B mesons and we believe that for these decays a lattice calculation of Fy tJ 4 (x 7 ) 
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FIG. 12: Point-like (pt), structure-dependent (SD) and interference (INT) contributions to the 
decay 7 r —>• £ 1 / 7 . The first (second) row corresponds to l = e (£ = /i). 


for a range of values of will prove to be very useful as a check of the range of validity 
of the point-like approximation. As stressed in the main body of the paper, such a lattice 
calculation, starting from Euclidean correlators is indeed possible. A new feature in the case 
of 5-decays in particular, one which is a consequence of the heavy-quark symmetry, is that 
the 5* and B are almost degenerate (m* B — tub — 45MeV). The radiation of a relatively 
soft photon can therefore cause the transition from a 5-meson to an internal 5* close to its 
mass-shell. Lattice calculations of the form factors would allow us to investigate the effect 
this small hyperfine splitting has on the size of the structure dependent terms as a function 
of A E. 

In the absence of lattice calculations of the form factors, we note the phenomenological 
analysis of Ref. (37|, based on the extreme assumption of the single pole dominance, B* for 
Fy and Bi{ 5721) for Fa (in reality many other virtual states contribute to the form factors): 

C V r, . C A 


5y(x 7 ) ~ 


Fa(x^) ~ 


1 + rn B 1 (5721)/ m B ’ 


(B17) 


x 1 — 1 + m 2 B */m 2 B ’ ' x 1 

with Cy = 0.24 and C A = 0.20. The corresponding ratios R± are shown in Figure fl4l from 
which it can be seen that under this assumption the structure-dependent contributions to 
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/?,[K-» ev(y)] 
0.00 I - ■ - ■ - 



AE(MeV) 


[K-> fiv(y)] 




FIG. 13: Point-like (pt), structure-dependent (SD) and interference (INT) contributions to the 
decay K —>• £wy. The hrst (second) row corresponds to l = e (l = /i). 


50 100 150 200 250 


50 " 100 150 200 250 


-1000 _1500_ 2000 2500 


FIG. 14: Structure-dependent (SD) and interference (INT) contributions to R\ for the decays 
B —> iwy. Going from left to right, the plots correspond to £ = e, i = fi and £ = r respectively. 


B —> ez / e 7 for E 1 ~ 20 MeV can be very large, but are small for B —> iw^f and B —> tu t 7 . 
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